Acoustic emission before avalanches in granular media. 
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Avalanches of granular media are mainly characterized by the observation and the measurement of 
the main angles of avalanche corresponding first to the movement of a group beads and to the whole 
movement of a great part of the grains. These characteristics do not give any information about the 
rearrangements of the grains inside the granular beads layer. As any movement of a grain produces 
a deformation of the structure, it is quite normal to expect a sound that will propagate inside the 
granular medium. We present an acoustical experimental study of the of avalanches precursors 
signals in spherical granular beads glass and silica aerogels powder (size of grains less than 80 
micrometers). Movements inside the material produce sound pulses, these acoustic emissions have 
been recorded with two piezoelectric transducers placed on the lower part of the material layer. Our 
results show clearly signals before any movement on the upper part of the bead's layer: recursive 
precursors in glass beads and non-recursive precursors in aerogels, for angles less than the first 
avalanche angle. Our methodology allow us to obtaining information of all internal movement in 
the granular material including avalanches, precursors, and small arrangements. 

PACS numbers: 45.70.Ht,43.60.+d 
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INTRODUCTION TO THE ACOUSTIC 
PROPAGATION IN AVALANCHES 



The great interest of the study of avalanches is oriented 
to the prevention of natural disasters and geophysical 
problems. Experimental works made by [l|, Q define an 
avalanche as a self organized critical phenomenon. When 
an avalanche is triggered by the inclination of the pile up 
to the threshold of instability, the free surface of the pile 
defines a critical angle, frequently called maximum sta- 
bihty angle Om- During the shde, a decrease of the slope 
of the free surface is observed until the avalanche stops. 
At this point the free surface defines a second critical 
angle: The repose angle Or- While the packing is being 
tilted, one can observe many superficial rearrangements, 
for angles just below 6m- So an avalanche occurs when 
the critical angle of the bulk of material 9m is passed 
by his own tilt and the pile "trying to organize itself 
slips material to keep the angle of repose Br- Recent 
works show evidence of events that occur before large 
avalanches [l|, [3, H, @, these events could be interpreted 
as avalanche predictors, they occur with angular regular- 
ity during an avalanche process. This angular regularity 
opens the prediction possibility of larger intensity events 
like larger avalanches. 

In [l| Nerone et al. studied the phenomenon of 
avalanches through an optical system recording a video 
of glass beads movements enclosed in a box when they 
change the angle of tilt of his assembly. The system 
counts the number of spheres that move when they 
change the angle of the box in one axis, achieving and 
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identifying precursors. This work defines the events that 
occur in a avalanche as: Small arrangements: small clus- 
ters formed on the free surface areas of the bulk when the 
small arrangements start, some of these beads are moved 
to other near place, and always this relation destabilizes 
the others beads. These events occur throughout all the 
experiment and anywhere on the surface. Precursors: 
These events involve all beads on the surface. It is noted 
that the beads are uniform distributed in the free surface 
because the precursors are not localized events. In this 
case, all beads can be simultaneously moving mostly in 
one step. They occur for angles 6 > 17° and have a reg- 
ularity of A9 ~ 1.8°. Avalanches: Events that begin as 
small perturbations located anywhere in the free surface 
with rapid growth and that will be propagated in some 
lower layers. Avalanches start at 9 > 27°. Q Study the 
precursors by a non-linear acoustic probing methodology 
which is sensitive to the state of weakest inter-grain con- 
tacts. For his experiment they emit a signal by an acous- 
tical transducer. This signal travels inside of the beads 
bulk in the movement box process in order to obtaining 
in other transducers the result signals with the informa- 
tion of all mixed systems. The measures are numerically 
handled in order to separate the bead's movement sig- 
nals and the noise, obtaining by this methodology the 
arrangement signals. This technique provides informa- 
tion of the state of weak contacts and only the small 
arrangements. 

We propose a different system to those used it by [l|, Q 
by an acoustical methodology to detect all the informa- 
tion of the beads movement, in a more simple way. An 
acoustical study of the pre-avalanche dynamics in a 3D 
granular packing in two different granular media enclosed 
in a box. We focus in the behavior of events that precede 
the avalanches and a non noise system. Any mechanical 
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variation or disarrangement of the initial state imposes 
constraints; these constraints have a vibrating signature 
that can be detected as an acoustical signal and can be 
recorded. The principal idea of our set up is to obtain 
information of the internal beads movement and not only 
the superficial movements. For this, we design a mechan- 
ical system formed by a transparent plastic box contain- 
ing the granular material. In the bottom part of the box 
2 transducers make the signals registration. For trigger- 
ing the avalanches the box changes softly the angle of 
tilt in order to register the signals produced by the ma- 
terial movement. We design experiments in two differ- 
ent granular materials in order to compare the avalanche 
phenomenon. 

The first experiment was carried out with glass beads 
of 3mm of diameter. When we varying the inclination 
angle, beads start to move and this movement will pro- 
duce a signal product of contact and friction between the 
beads of glass that move inside the material. The first 
movements in this assembly may be seen as precursors of 
the avalanche studied by [3, S S 01 ■ These precursors 
will be produced by structural instability's inside of the 
arrangement of the beads by the inclination of it. 

The second experiment was performed with silica aero- 
gel. Snow is very difficult to manage, for temperature 
problems and non reproducibility factors; the question 
has been to identify a material with similar snow char- 
acteristic, especially with similar cohesion properties. It 
was already demonstrated by Ishida [9] that snow can be 
treated explicitly as a granular-porous material. Previ- 
ous studies on granular silica aerogel have provided 
extensive information on the characteristics in very low 
density and highly porosity materials. When the size of 
beads is less than 40 micrometers in granular materials, 
[l^ gravity force is less intense than the electrostatic. As 
a result, in aerogel powders (diameter size between 40- 
80 micrometers) , electrostatic and gravity forces compete 
each other. Following this property, the cohesion forces of 
snow can be replaced by the electrostatic forces in aero- 
gels, emulating the snow behavior. This powder gives the 
same "natural" features; fractures and block avalanches 
producing a perfectly mimicking like in a snow at a very 
small scale. The speed of sound is very low in silica aero- 
gel powder (around 60 m/s in comparison to 250 m/s for 
snow) ll] however the same feature is also expected for 
the fractures propagation acoustic emission signals. In 
both cases we make the signals frequency analysis. This 
analysis provides a way to view the avalanches process 
as a statistical phenomenon that gives information of the 
statistical periodicity during the avalanches process. 

2. EXPERIMENTAL SETUP 

Our study tries to focus in conditions for measures 
without external noise in order to obtain clear acoustic 
signals. The assembly took place in a transparent plas- 
tic box of Methacrylate of 4cm Height, 40cm of long and 




FIG. 1: Scheme for the miniature avalanches experiment. The 
angle of the box 9 will vary softly when the ball is deflated. 
The ball touch the box through a metal rod in order to de- 
flate the ball avoiding mechanical noise. The velocity of the 
downing process is controlled by the air valve 



40cm wide (see FIG 2). The height of the layer was 3cm 
for glass beads and 2cm for the silica aerogel powder. 
The box changes its tilt angle using a mechanical system 
articulated on its " upper" side, while the lower side is mo- 
bile and the whole box can rotate around this side axis. 
The rotation is controlled through an inflated balloon. 
When the balloon is slowly deflated the box is rotating 
slowly and its position and velocity can be manipulated 
by stopping or regulating the quantity of air that's let 
out in the deflating process with a valve. The pile is 
softly tilted to avoid noise disturbance by the mechani- 
cal system, up to the threshold of instability where the 
avalanches are triggered. On the box bottom plate two 
Piezo-electric transducers of resonant cells (resonance fre- 
quency around 3 kHz) make the acquisition. The acous- 
tical signal has been recorded on both transducers at a 
12 bits quantization oscilloscope for a duration of 50sec, 
400Hz sampling frequency and 20k samples for the sil- 
ica powder, and 200sec, IkHz of sampling frequency and 
200k samples in glass beads. We can see a Description of 
the experimental set-up in FIG (1): 

In the lowering process the velocity of angle variation 
is roughly linear because the angular velocity is slow. 
The angle variation can be directly related to the box 
lowering time in order to interpret the angular variation 
in the graphics. 

This simple system allows repeatable experiments with 
the advantage of minimal noise. In order to minimize 
the mechanical wave conduction between the set up and 
the laboratory floor, we place the experimental set up 
over a foam surface. FIG (2) represent the box rotation 
without any material inside as a test of noise in the lab. 
The signal level was around 0.00079mV while the small- 
est recorded signals corresponding to an avalanche have 
amplitude of around 0.02781mV and the regular signals 
have amplitude between 0.5mV. In the FFT graphics for 
frequencies domain, we report two specific frequencies in 
50Hz and 232 Hz. These frequencies are present in all the 
experiments and correspond to the lab noise produced by 
a motor noise in the nearest of laboratory building. 



3 




16 



30 



FIG. 2: System downing measures without material. Middle Transducer 500Hz 10k samples 20sec. And FFT 
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FIG. 3: Beads of glass 3mm. Middle Transducer IKHz 200k samples 200sec and FFT 



We perform 10 experiments for glass beads and 10 ex- 
periments for aerogel powder. Each one under the same 
{'xpaiincntal conditions. Wc compute the angular veloc- 
ity in the box downing process for each experiment using 
the acquisition time and the familiar expression of lineal 
angular velocity v = 0/t in order to obtain the angular 
path for each case. The experimental time was assumed 
long for be able to approximate at a lineal angular ve- 
locity behavior. With this velocity wc can compute the 
signals angular position for each graphic. 

In the FIG (3 and 4) we present the signals obtained 

for the transducer located in the middle and upper part 
of the Plexiglas bottom box in one of the 10 experiments. 



The signals observed correspond to the beads movement 
inside the box. The amplitude of these signals is associate 
to the quantity of moving beads. Less signals events are 
associate to the precursors and more amplitude to the 
larger avalanches. There are a remarked different be- 
tween both kinds signals in the graphics. Clear precur- 
sors events signals are observed with a angular regular in- 
terval mean of = 2.09°. The compute of the standard 
deviation (SD) of these mean value give us a measure of 
the data dispersion and the proximity in the regularity for 
these events. More similar events will be produce small 
values in DE. In this experiment ED = 0.042. The pre- 
cursors start for an angle of ^ = 15.8° and the avalanches 
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FIG. 4: Beads of glass 3mm. Upper Transducer IKHz 200k samples 200sec and FFT 
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FIG. 5: Beads of glass 3mm. Upper Transducer IKHz, 200k samples 200sec. Precursor Start Point, Larger avalanches and A9 
value for 10 experiments. 



start in 6 = 27.53°. For all the 10 experiment we present 
the table 1 and FIG (5). In these we can observe the 
behavior in 10 ten experiments in order to compare. We 
present the measures for each experiment for: the precur- 
sor's start angle, mean precursors interval, mean precur- 
sors interval SD, and avalanches start angle. We present 
the arithmetic mean for the characteristic angles for all 
experiment (see table I). These values present a clear ten- 
dencies in the precursors start angle 9 — 17.04° with a 
dispersion of SD = 0.326°, avalanches start 9 = 28.41° 
with SD = 1.29°, the size of the precursors mean in- 
terval was A0 = 1.94° with SD = 0.18°. These values 
are agree with the Nerones work in [l| where they values 
for the precursors start in 9 > 17° large avalanches for 
9 = 27.3° and precursor angular interval of A9 — 1.8°. 



The small differences are maybe associate to the fact that 
the avalanches are not a lineal phenomenon, for that its 
impressive to make a more higher number of experiment 
in order to compare the mean behavior of it. 

The Fourier transform gives us information about the 
events repetition inside the signal. If we imagine one 
beads downing in a beads layer. This bead can colli- 
sion against the others and produce a sound that can be 
register in a signal. These signals will have a period as- 
sociate to the time between collisions, and this period a 
frequency associate that can be observed in a FFT. For 
our case we have a lot of bead moving in the box, and 
the FFT spectrum gives us information of the statistical 
repetitive events. We present in the FIG(6), 4 experi- 
ments with the same conditions, a clear Gaussian behav- 
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TABLE L 3mm Glass Beads information for each experiment. 




FIG. 6: 4 Glass beads experiments EFT. Middle transducer 
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TABLE IL Aerogel information for each experiment. Middle transducer. 
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TABLE III: Aerogel information for each experiment. Upper transducer. 
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FIG. 7: Aerogel. Middle Transducer. 400Hz 20k samples 50sec and FFT 
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FIG. 9: 4 aerogel experiments FFT. Upper transducer 



ior cantered in approximate lOOHz is reported. These 
frequencies have a period of O.OOlsec. We assume this 
period as a statistical measure of time collisions in the 
bulk of material when they change his tilt. These sig- 
nals correspond to those reported for Nerone in [l[ as 
small arrangements in his optical system. This acousti- 
cal signals methodology give a easier and complete way 
to obtain all the avalanche informations and to apply in 
field observations. 

In the aerogel the results are different. The avalanches 
are block avalanches like in the snow caseQ- This 
avalanches start by internal fractures inside the material, 
our propose is to detect the fracture signals as a precur- 
sors events. Is converdent to remember that the aerogel is 
a damping material [lO| . This feature prevents the detec- 
tion of signals coming from distant source locations. As 
a result, there are differences recorded amplitude values 
between the transducer located at the top and the mid- 
dle. In our case we perform 10 experiments. Each one 
with similar characteristics, in a 50 seconds measurement 
interval 20k samples and 400IIz. FIG. (7 and 8) show the 
behavior of an amplitude and frequency of one of the 10 
experiments in aerogel in the middle transducer and up- 
per respectively. The amplitude of the signals in Figure 
7 are different from those show it in FIG (8) for the rea- 
sons explained above, however, the Fourier spectrum in 
both cases give similar results 167IIz and 74Hz signals. 
These frequencies correspond to short periods vibration 
caused by the attenuating properties of the material. In 
the FIG (9) we present the behavior of the frequencies 
in 4 different experiments. The 50IIz frequency is not 
taking in a count in the analysis because correspond to 
the laboratory noise as we can clearly see in the graphics 



in FIG 3. 

As a result of the graphics observation we can affirm 
that does not exist angular regularity in the signals am- 
plitude events in FIG (7 and 8). This is reaffirmed in 
the table II and III where the results present a higher 
SD values in comparison of his own mean for the pre- 
cursors interval A0„ for each particular experiment. In 
many cases there are not presence of significative ampli- 
tude signals for consider as avalanches or even precursors. 
However this situation only occur in 1 of the 2 transducer 
(see tables II and III). If there are not signals in one, the 
other yes have it. That's because there are not enough 
energy to propagate the acoustic wave for all the material 
bulk. 



3. CONCLUSIONS 

Our low noise levels system allowed us to obtain the 
signals that occur when the pile is tilted to the threshold 
of stability in beads of glass and powder of aerogel. In 
the glass beads case, we reproduce the results obtained in 
[l| using a acoustical system and comparing in 10 exper- 
iments. We processing the signals with a FFT obtaining 
the frequency spectrum. The frequencies analysis show 
us results associated to the movement of the complete 
bulk of material. Obtaining by this methodology the 
mean of the time movement period in the glass beads or 
an mean of frequency movement in small arrangements. 
In the aerogel powder, it's the first evidence of block 
avalanches precursors signals. In this case we obtain 
signals only in the proximity of transducer because his 
damping characteristic. However middle intensity sig- 
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nals arc reported before the great avalanches, showing 
evidence of non regular precursor. This system confirm 
the possibilities of use this kinds of acoustical setup in 
real field scenarios. This signals despite not have re- 
cursively can be take it as a predecessor of the larger 
avalanche events. Despite both are a granular material 
the behavior of avalanche granular dynamics are differ- 
ent. The interaction between the grains produce different 
phenomenon because the different individual weight for 
the glass beads particles and the aerogel particles. In the 
glass beads case (heavier material) the gravity produce 
over the grains a faster individualized movement, the 
gravity force is more significant in comparison of cohesion 
forces. This force act over each particles un-stabilizing 
the material bulk when it arise his critical angle. This 
instability produce the localized phenomenon of precur- 
sor and the larger avalanche. The precursors are caused 
by the dis-stabilisation of the individual grains over his 



neighbor. When the instability arise a more significant 
level, all the material slip down. In the aerogel this in- 
dividual behavior is more retarded because the electro- 
static force interaction inter-grain is more significant in 
comparison to the gravity and this act like the cohesion 
force. The grain contact keep without moving the aerogel 
bulk particles even near to critical angle, the avalanches 
occur only when the threshold of instability is more in- 
tense. The avalanche start in fractures, when the first 
particles are moved because the effect of gravity force 
wining vs the cohesion force, the material start to sep- 
arate this cause a chain reaction that separate a lot of 
near particles causing and propagating the fracture. If 
this instability continue the fractured material is divided 
in small piece of materials, and the gravity slip down 
this more heaviest pieces increasing the instability and 
producing the avalanches. 
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